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Abstract 

In this paper we consider the problem of determining the law of binary stochastic processes 
from transition kernels depending on the whole past. These kernels are linear in the past 
values of the process. They are allowed to assume values close to both and 1, preventing 
the application of usual results on uniqueness. More precisely we give sufficient conditions for 
uniqueness and non-uniqueness. In the former case a perfect simulation algorithm is also given. 

keywords: Perfect simulation, Coupling, Chains with complete connections. 
AMS classification: 60G99, 68U20, 60J10. 

1 Introduction and main definitions 

In this paper we consider processes indexed by Z taking values in a finite alphabet G. These 
processes are constructed from a transition kernel which depends on the whole past, i.e. a map 
p : G x G~ N+ — > [0, 1] such that, for any choice of wZ^ = (u>-i, W-2, ■ ■ •) S G~ N+ , p(-|wl^ ) is a 
probability measure on G. For the basic setting of the theory of these processes we address the 
reader to |CFF02| where an historical overview is also given. In particular we recall the following 

Definition 1. A G-valued process X = {X n ,n G Z} defined on an arbitrary probability space is 
compatible with the transition kernel p if for any g £ G 

P(X n = g\X^) = pGHX!£), (1) 

almost surely, for any n€Z. 

It is clear that the compatibility condition refers only to the law of X, so we are allowed to 
speak about compatible laws, whose set will be denoted by Q{p). When \G{p)\ = 1 the law of any 
compatible process X is always the same in which case we say that uniqueness holds for p. 

It is immediately seen that G(p) is convex. Moreover, if the kernel p is continuous w.r.t. the 
product topology in the second coordinate, then G(p) is non empty (see |FM05| ). In this paper we 
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will establish sufficient conditions for uniqueness for a certain class of continuous kernels, giving 
in these cases an explicit construction for the unique compatible process. 

In the sequel we will consider binary processes taking values in G = {— 1, +1}. The transition 
kernels we are interested in have the form 

_. oo 

/G?|w:L) = 2 + fE^-*> sec, wI^GG-^, (2) 

k=l 

where 9 = {9 k, k E N+} is a sequence of constants with the property X^fcLi l^fel = 1- A kernel of this 
form is clearly continuous, since by specifying some large but finite portion of wZ^, the oscillation 
of p s (g\-) can be made arbitrarily small. Moreover, since p 9 (g\wZ^ c ) = p 9 (— g\ — wzjo), Q{p e ) 
is invariant under the transformation changing the sign to the whole configuration. Therefore, if 
= 1 an d {X n } is the unique compatible process, P(X n = ±1) = 1/2 for any integer n E Z. 

Binary linear kernels were first introduced in [CFF02J, as simple examples of transition kernels 
depending on the whole past. More precisely they added a "bias", that is an extra term with 
#o £ [ — 1, 1] ; if a bias term is present the model is well defined provided YlT=o \@k\ < 1- But then 
they got YlkLi \@k\ < 1) i n which case uniqueness always holds. In fact the unique compatible 
process can be constructed by means of a simple perfect simulation algorithm, to which we will 
briefly refer in the last section. Alternatively, as in |FM05| . in this case one can prove uniqueness 
by means of Dobrushin's contraction condition. 

Thus only when 9q = the choice X^fcLi l^fcl = 1 ^ s allowed. As a consequence the convex hull 
of the range space of p s (l\-) (and of p s (— 1|-)) is the whole interval [0,1]. Therefore there is no 
uniform positivity condition for this kind of kernels, preventing the applicability of general criteria 
for uniqueness (sec [FM05]). Indeed uniqueness does not always hold. The simplest situation in 
which uniqueness does not hold is when 9k > for all k E N+. Then it is easy to check that 
the two constant processes constantly equal to +1 and —1 are compatible. More generally we can 
establish the following 

Theorem 1. Suppose that 

a) either 9] t >Q for all k E N+, 

b) or 9^ < for k odd and 9^ > for k even. 
Then uniqueness does not hold for p e . 

On the positive side, our main result is the following sufficient condition for uniqueness, which 
will be supplemented by a perfect simulation algorithm producing the unique compatible process. 
In order to state it, we set Ag = {k E N+ : 9^ ^ 0}. 

Theorem 2. Suppose that all the following assumptions hold: 

i) gcd(A e ) = I; 

ii) there exists either an even integer m E N+ with 9 m < 0, or an odd integer mi with 9 mi < 
and another odd integer mi with 9 m , 2 > 0; 
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Hi) Efcll(Em=fc \ e m\) 2 < +00. 

Then uniqueness holds for p e . 

We briefly comment the relation among Theorem Q] and Theorem [2J Condition i) in Theorem [2] 
is necessary for uniqueness (Proposition [1]); condition ii) is the negation of assumptions a) and 
b) of Theorem [TJ At the present time only condition iii) is preventing the statement of a set of 
necessary and sufficient conditions for uniqueness. We will briefly comment about this within the 
conclusions. 

In Section 2 we will prove Theorem Q] and related results. In Section 3 we will prove Theorem [2] 
and related results. In Section 4 we give the simulation algorithm explicitly with a few conclusive 
remarks. 

2 Coherent sequences and nonuniqueness 

The set A = Ag allows to define an oriented graph structure T = (Z, A) over the set of integers. 
There is an arc (n, m) G A if m < n and n — m G A. Such an arc is marked by the sign of n -m- 
It is clear that, by its very definition, an arc in T can be always shifted, since if (n, m) G A, then 
(n + k,m + k) £ A, for any k G Z. An oriented path is a finite sequence {(nj,nj + i) G A, i = 
0, — 1}. This path joins no with n;. Each oriented path has a sign, which is given by the 
product n!=o s ig n ($ni-n l+ i) °f the signs of all its arrows. 

By neglecting the orientation of the arrows in A we obtain the set of edges A of the undirected 
graph r which is called the symmetrization of T. In this framework a path is a finite sequence 
{{ni,n i+1 } G A,i = 1, ...,£}. The sign of the path is given by n!=o sign(0| n ._ n . +1 |). Here the 
endpoints no and ni joined by the path could coincide, but in this case we prefer to use the 
word cycle rather than path. By convention for the path we always take no > n;, unless when 
specified otherwise, but we say that an edge {nj,nj + i} is traveled in the opposite direction when 
(nj + i, n,j) G A, therefore rij < Wi+i. By ordering the edges of the path {{n^, nj + i} G A, i = 1, . . . , 1} 
by decreasing rather than increasing the index i the same path is traveled in the direction which 
is opposite to the previous one. Finally, as arcs can be shifted in A, a whole path can be shifted 
in A in a natural way. 

By Bezout's identity the group generated by A is gcd(A) x Z, and the connected components 
of r coincide with the congruence classes mod gcd (A). In particular, when gcd(^4) = 1 the graph 
r is connected. This means that for any pair of sites n, m G Z there exists a path joining them. 

Proposition 1. The condition gcd(vl6i) = 1 is necessary for \G(p e )\ = 1. 

Proof. Suppose gcd(vl) = ko > 1. Define 9* k = 9kk for k G N+. Let X* be any p e * -compatible 
process. It is easily seen that any process X = {X n ,n G Z} such that X^ 1 ) = {X n ) :o+ f l ,n G Z} 
has the same law of X*, for any h = 0, 1, . . . , ko — 1, is a ^-compatible process. Since the joint 
distribution of X^ for h = 0, 1, . . . , ko — 1 is arbitrary it is clear that uniqueness does not hold 
for p s . In particular one can take all X^'s either equal or independent, in order to produce two 
different elements of G(p 8 )- Q 
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Another source of nonuniqueness comes with the following definition. 

Definition 2. A configuration s E G z is coherent with the sequence of coefficients 9 G G N+ if for 
any (n, m) E Aq it is s n s m 9 n -m > 0. 

Coherent configurations come always in pairs, since if s G G z is coherent, the same is true for 

— s. 

Proposition 2. A configuration s E G z coherent with 6 is a degenerate compatible process (i.e. 
$s G Q{p 8 ))- Thus, when a coherent configuration exists, then \G{p d )\ > 1- 

Proof. First notice that if s n s m Q n - m > 0, then s n s m 6 n - m = |0 n _ m |. Now, if X„ = s n , n E Z, 
then, for any n E Z, it is obviously P(X n = s n |X™"^ = s™"^) = 1. But on the other hand 

1 = ^S^ = ^iE^ ( 3 ) 

k=l k=l 

which agrees with ([5]). □ 

The following result characterizes the existence of coherent configurations in a convenient way. 

Lemma 1. If a coherent configuration exists all the paths joining any two given sites have the 
same sign. Conversely, if all the oriented paths joining any two given sites have the same sign, a 
coherent configuration necessarily exists. 

Proof. First notice that, if a coherent configuration exists, for any path {{n,,nj+i} E -4,}i=o,...,i-i 
we have s m s ni n!=o @\ni-n i+1 \ > 0> from which necessarily all other paths joining no and n/ must 
have the same sign. Conversely, choose a representative for each connected component of T, and 
define arbitrarily the configuration s over these sites. Suppose no is one of these representatives 
and consider any path {{nj,nj + i} G A}i=o,...,i-i- By defining 

11 

s ni = s no JJsign(0| ni _ n . +1 |), (4) 

i=0 

one can extend unambiguously s to a coherent configuration. It remains to prove that if all the 
oriented paths between any two sites have the same sign, the same is true for all the paths. In fact, 
suppose that there are two paths with different signs joining two distinct sites. By following one 
of the two paths in the opposite direction and concatenating it with the first we get a cycle. By 
assumption the sign of this cycle is negative. We can now produce two oriented paths joining the 
same endpoints by shifting and concatenating the arcs which are traveled in the same direction. 
These two oriented paths necessarily have opposite signs. □ 

By the previous result, when gcd(A) = 1 and all the paths between any two sites have the same 
sign, a coherent configuration s is uniquely determined by the choice of the sign sq at the origin. 
Next observe that by shifting a coherent configuration, another coherent configuration is obtained. 
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Now if a configuration is invariant under the shift, it is necessarily constant. If it changes sign under 
the shift, it is necessarily one of the two "checkerboard" configurations alternating periodically — 1 
with +1. Thus we have the following result. 

Proposition 3. When gcd(Ae) = 1 a configuration coherent with p e exists only in the following 
two cases: 

a) either 9k > 0, for all positive integers k; 

b) or 9k < for k odd and 9k > for k even. 

Proof. The two cases correspond to the two possible pairs of coherent configurations, the constants 
and the checkerboards. It is easily checked that the constant configurations are coherent in case 
a), and the checkerboards in case b). □ 

From Propositions Q] [2] and [3l Theorem [1] immediately holds. 

Before turning to discuss sufficient conditions for uniqueness we observe that any pair of con- 
figurations ±s E G z is equivalently described by an unordered partition of Z into two subsets, 
which are the inverse images of the two possible values —1 and +1. 

If gcd(.A) = 1 and the checkerboard configurations are coherent, this partition has the property 
that each pair of sites with the same parity is joined by a positive path, and each pair of sites with 
different parity is joined by a negative path. If gcd(A) = 1 and the constant configurations are 
coherent, then we have a degenerate partition with a single class: all sites are joined by a positive 
path. This interpretation will be useful when dealing with partitions of this kind for suitable 
subgraphs of T. 

3 Random subtrees and uniqueness 

We start by observing that sampling from p s given in ((2J) is realized by drawing a random variable 
K with P(K = k) = \9k\-, for k = 1, 2, . . . and returning X* = sign(#/^ )u>_^. Indeed, by setting 
i/k = sign(9k)w-k, one has 

1 1 1 00 

p(x* = i|w^j= £ \e k \ = -{ N + C 1 - E N)} = 2 + 2^ fcl?/fe ' (5) 

k:y k =+l k:y k =+l k:y k = -l k=l 

and since \9k\yk = QkW-k, the statement is proved. 

Thus, given a p -compatible process X = {X n ,n £ Z} we have the following resampling scheme, 
that does not change its distribution. Given any site n € Z, we select K n from the distribution 
{9k, = 1,2,...}, and then pick the arc (n, n — K n ) from the set of arcs in A outcoming from n. The 
resampled value X* = s\gn.{9K n )X n _K n is then obtained by sending the value at the sink n — K n 
of the chosen arc back to the source n, changing the sign if the sign of the arc is negative. 



5 



By drawing K n as above, independently for each n £ Z, we obtain in this way a random 
subgraph r r = (Z, 0) of T = (Z,^4), and we call its symmetrization T r = (Z,23). This will play a 
fundamental role for uniqueness. Since in B there is only one outcoming arc from each site rGZ, 
in B there is no cycle: thus the connected components of T r are trees, and such a graph is a forest. 

Let us define 

Jf = a {Ki :k<i<m), (6) 

where — oo < k < m > +00. Define the left tail u-algebra F-oo = ^m&F^oo- The number N of 
trees which forms T r is a random variable is .F-oo-measurable since it is possible to read N from 
the values {K n ,n < m} no matter how far the site m € Z lies in the past. In fact all these trees 
necessarily "have to be born at —00", since each oriented path in B can be extended infinitely to 
its left. Since the K n J s are i.i.d., it follows that N is a.s. constant. In the following we will prove 
that provided the assumptions of Theorem [2] do not hold (otherwise coherent sequences exist) and 
N = 1 a.s. (which clearly implies gcd(vl) = 1), there is uniqueness for 9. 

Before proving this result we exhibit what in these cases will turn out to be the unique com- 
patible process. Since N = 1 any two sites will have a unique path connecting them. Hence we 
can partition the graph T r in two disjoint components in the following way. Starting from an 
arbitrary site, say the origin, we define the random set of sites S which are joined to the origin by 
a positive path (including the origin itself). Each pair of sites belonging to this set is joined by a 
unique positive path. The same property holds for the complement S c . Moreover, each site in S 
is connected to any site in S c by a negative path. In order to define a process from this partition 
it remains to select a random sign common to sites in S (the opposite one is given to those in S°). 
This is achieved by assigning to the origin one of the two signs with equal probability. To sum up 
we define 

X n = Y(2-l s (n)-l),neZ (7) 
where Y is independent of S, with distribution P(Y = ±1) = 1/2. 

Proposition 4. The process X = {X n ,n G Z} defined in ([?[), is compatible with p 9 . 

Proof. It is clear that, since the origin is connected to any site almost surely, also the sign of 
such a site will have the same symmetric distribution. Since the unordered partition {S, S c } is 
invariant by translation of the K^s, the process {X n ,n 6 Z} is certainly stationary. Thus, to 
verify compatibility, it is sufficient to prove that resampling the value in 1 according to ([2]) will 
not change its distribution. For this we substitute Kq with Kq having the same distribution, 
independently of all the K n 's. This means that the subtree attached to the origin, which originally 
is attached to the site — Kq, becomes attached to the site —Kq. The resulting partition is the 
one which would be obtained directly from {K n , n 7^ 0, Kq}, hence it has the same distribution as 
before. Moreover the sign of the origin is either that of X-k* } or the opposite, depending on the 
sign of Ok*. But the sign of X_k* has again a symmetric distribution, independently of all the 
{K n ,n 7^ 0, Kq}, because of the way Xq = Y was defined. This ends the proof. □ 

For proving our uniqueness result it is now necessary to introduce the processes X _n,w , 
parameterized by n £ Z and w G G~ N +, with the "boundary conditions" = w-i for 



6 



i = 1,2,... and for all the remaining sites m > —n obtained by following the unique oriented path 
{(Ni, Ni+\),i = 0, . . . , L m — 1} starting from iVo = m in T r (which means that N^i = N — Kn { ), 
stopped when it "reaches a boundary site" (which means that L m = inf{i : Ni < — n}). Finally 



so that X _n ' w is obtained from the boundary conditions by following iteratively the sampling 
scheme introduced at the beginning of the section. All these processes are thus constructed on the 
probability space where the process {K n ,n G Z} lives. 
Here is the promised result. 

Theorem 3. Consider a transition kernel p of the form p|) with the following assumptions on 
the coefficients: 

a) there exists either an even integer m G N+ with 9 m < 0, or an odd integer mi with 9 m < 
and another odd integer mi with 6 n > 0; 

b) the random graph T r is a.s. connected. 

In order to prove the above result, we prove first a very useful lemma. 
Lemma 2. Under conditions a) and b) 



Proof. First observe that assumption a) ensures that there are no coherent configurations. By 
Lemma Q] we know that this requires the existence of two oriented paths in T with the same 
endpoints, say and —niQ < 0, but with different sign. By excluding a common initial subpath 
to these paths, one can reduce tjiq to the minimum possible value. Next denote the positive path 
by 7 + and the negative one by 7~. We write p+ (resp. p_) for the probability that 7 + (resp. 7~) 
occurs, an event which is J^^-measurable. 

Let us introduce the random times {T,,i G N + }, defined by recursion from To = 0. These 
random times are a subset of the endpoints of the arcs of the infinite oriented path {(iVj, A^+i), i G 
N} in r r , starting from iVo = 0. Whereas Ni — Ni+i = K^, for % G N, the Tj's need to have a 
distance not smaller than mo, in order to give to 7 + and j~ a positive probability to happen. This 
explains the following definition. 




(8) 



i=0 




(9) 



Tj + i = sup{m < Tj — mo : there is a path in T r joining Tj with m}. 



(10) 



Furthermore, for i G N+, define the marks 
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Finally define J n = sup{m : T m > — n}, = \{i < J n : Y, L = +1}| and R n = \{i < J n : Yi = — 
It is clear that lim n ^oo(i?~ + i?+) = oo almost surely. Therefore there exists M n ~\ oo and e n \. 
such that P(R~ + R+ > M n ) > 1 - e n . Then 

oo 

P(X~ n ^ = +l)< £n + P(X n ' w = +l\Rn +R+ = k)P{R- + R+ = k) 

k=M n 

<e n + max P(X(7 n ' w = +1\R~ + R+ = k). 

k>M n 

Since, each time that a segment which is a shift of 7 + is replaced by the same shift of 7" the value 
of XJ~ n ' w changes sign we can bound the previous formula by 

< e n + max \ max.{P(R~ is odd|i?~ + R£ = k),P(R~ is evenji?" + R^ = k)} > . (12) 

k>M n I J 

In order to bound the previous expression we notice that R~ conditioned to the event {R~ + R£ = 
k} has a binomial distribution, namely B(k, p ^ p ) • For any random variable Xj~ ~ B(k,p), with 
p € (0, 1), it is easily proved that P{Xk is even) goes to 1/2 as k — > 00. As a consequence we 
obtain that the expression (|12|) is bounded from above by 

1 

- + e n + d Mn 

where lim^oo 5k = 0. Since the same bound can be produced for P(X n,w = —1), this ends the 
proof of the lemma. □ 

Notice that, since for random variables X and Y taking values in {— 1,+1} with laws C{X) 
and C(Y) it holds 

\\C(X) - C{Y)\\ TV = sup \E(f(X))-E(f(Y))\ = \P(X = l)-P(Y = l)l 

\\f\\°°<i 

we can reformulate the above lemma as 



lim sup ' ) " £(X ' hW2 )\\TV = 0. 

l ^ 0C w i eG- I '+,*=l,2 



Proof of Theorem In order to prove uniqueness let us consider two p e -compatible processes X 
and Y defined on the same probability space. By using a sequence {Ki,i € Z} of i.i.d. random 
variables with distribution 8, independent of X and Y, it is possible to sample recursively all 
the values at sites m > —n, according to the sampling scheme described at the beginning of the 
section, getting process X* and Y* that have the same law of X and Y, and that will be identified 
with them from now on. Notice that the invariance of the laws under resampling allows to write 

„ v— n — 1 „ -\r — n — 1 

x = x -n,X_ oc )Y = X- n ' Y — . 
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Since N = 1 a.s., for any finite window of sites [I, L] = {i £ Z : I < i < L}, with — oo < / < 
L < oo, the root of the smallest subtree of T r containing [I, L] is a proper random variable T^^y 
For any function / : Q L ~ l + l — R ) with ||/||oo < 1 we write 

£{/(Xf ) - /(Yf)} = £{! {%i]> _« } (/(Xf) - /(Yf ))} + S{l {T[l£] <_n } (/(xf ) - /(Yf ))} 



s 



< £{1 {T[ , £]> _ 5 } (/(Xf) - /(Yf ))} + - (13) 

for any e > 0, for n sufficiently large, since Tn^ is finite a.s. 
Observe that is Fk L ,+i -measurable, and 



Xf - X T[l L] Zf , Yf - Y T[l L] Zf , 

where Zf is an (L — Z + l)-dimensional i]+1 -measurable random vector: moreover Xt^ l) and 
Yt [; l j are independent of J 7 ^ Therefore 

£{1 {T[! L]> _™ } (/(Xf ) - /(Yf ))} = E{l {T[lL]> _n }E (f(Xf ) - /(Yf)|^ (14) 

Finally, setting / z l(<?) = f(g x Zf ), for g £ G, we have 

£(/(Xf ) - /(Yf = E{/ Zf (X T[Iii] ) - / zf 0V,)}, (15) 

where E means expectation w.r.t. to the joint law of l] and L] . Thus, plugging the above 
expression into (fl4"j) . we get 



£{l {w _n } (/(Xf) - /(Yf ))} = E(l {T[i L]=t} K(f zt ,(X t ) - f Z L(Y t ))) 

-n/2<i<l 

< Y, nT { i,L]=i) sup \\C{X7^)-L{X7 n ^)\\ TV . (16) 

in^-^l Wi,W 2 
— n/2<i<l 

By Lemma [2] and the subsequent remark the total variation term can be made smaller than | as n 
is large enough. Plugging (|16p into (|13p we thus obtain that the total variation distance between 
the laws of Xf and Yf does not exceed e, proving the desired result. □ 

The natural question is under which conditions on the sequence of coefficients 9 i is possible to 
ensure that N = 1 a.s. or, equivalently, that any two sites, say e.g. and —1, are a.s. connected 
by a path in T r . We need to follow the two infinite oriented paths in T r starting from this two 
sites to check if they coalesce or not. This can be done by keeping track recursively in time of 
the distance between the two paths, using the rule that it is always the rightmost path which is 
updated with the addition of an arc. If the rightmost path reaches the site i, it jumps backwards 
to i — Ki, and so on, independently of the arcs chosen in the past. In this way the distance between 
the two paths follows the Markov chain 

Y n+1 = \Y n -K n+1 \,neN,Y = i (17) 
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where the 7Q,i = 1,2,... are i.i.d. with distribution P(K\ = j) = 9j. This process is known 
in the literature as the von Schelling process (see }Fel71j ). It evolves as a random walk with 
negative increments until it falls below zero, in which case it is reflected on the positive side. 
When gcd(Ag) = 1, this process is irreducible (see [Bou86j ) . Moreover we have the following 
result. 

Proposition 5. Consider the random graph {"L,B) with B = {{i,i — Ki},i E Z}, where the 
Ki,i = 1,2, .. . are i.i.d. with distribution P{K\ = j) = 8j. Assume gcd(/c : 9^ 7^ 0) = 1. Then 
the number of connected components N of (7*, B) is a.s. equal to 1 if and only if the corresponding 
von Schelling process is recurrent. 

Proof. It is clear that N = 1 a.s. if and only if the state is hit a.s. when the von Schilling process 
is started in 1. By irreducibility, this happens if and only if the process is recurrent. □ 

An invariant measure for the von Schelling process was written down in |Bou86| . Whenever 
it has finite total mass then the process is positive recurrent. This happens whenever E(Kq) = 
J2T=i P(K-o > k) < +00. Recently some null recurrent cases have been covered by a result 
in [PWllj . which is reported below. Incidentally, in the cited paper it is also given a class of 
distributions with very heavy tails for which the von Schelling process is transient. 

Theorem 4. IfY^kLi p ( K o > k) 2 < +00, then the von Schelling process is recurrent. 

By putting together the two Theorems [3] and 0] and Proposition [5] we have finally obtained the 
promised Theorem [2j 

We finally remark that N = 1 can be seen as a percolation type condition. For other examples 
of percolation techniques applied to perfect simulation in dimensions higher than one the reader 
is referred to |DSP08| . 

4 Perfect simulation and conclusions 

In the following we give a possible implementation of the simulation of the unique p -compatible 
process in a finite window [/,L].For the sake of simplicity and without any loss of generality we 
have set / = 1. 

Simulation algorithm for producing X = (X(l), . . . , X(L)) 

1. For i = 1, L, set X(i) = 1, S{i) = i, Y(i) = 1 and Z(i,j) = 5ij,j = 1,...,L. 

2. Draw K from the distribution 0. 

3. Compute J = argmax{S'(i)y (i)}. 

4. Set S{J) <- S(J) - K and X{J) <- X{J) x sign(^). 

5. If S{J) + S{i),i^ J go to 2. 

6. Else let / s.t. 5(7) = S(J), set Y(J) = and Z(J,I) = X(J) x X{I). 
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7- If Ef=i Y(i) > 1, go to 2. 

8. Else let H such that Y(H) = 1, and set X(i) = 0,i ^ H, drawing X(H) as a symmetric 
Bernoulli sign. 

9. Output X <(— Z^X, the product of the matrix Z L with the column vector X. 

A few words of explanation. The vector S = (5(1), . . . , S(L)) contains the positions of particles 
started from the sites of the interval [1, L], performing independent random walks until coalescence. 
These particles will be labeled with their starting sites. Their displacements are drawn from the 
distribution 9; at each step it is always the rightmost particle which moves to the left. The vector 
X = (A(l), . . . , X{L)) keeps track of the current sign of the paths followed by the particles. When 
particle J lands in a site already occupied by particle I we set to the J-th value of the auxiliary 
vector Y = (Y(l), . . . ,Y(L)). This means that particle J remains "freezed" in that site and does 
not play a further role, aside from the fact that in the auxiliary L x L matrix Z the (J, I)-th 
element registers if the paths followed by the two particles until coalescence have or not the same 
sign. This agreement or disagreement will be kept in the final vectors of signs. When the next to 
the last particle lands in the site Thn occupied by the last one, say particle H , by computing the 
L-th power of the matrix Z it is possible to reconstruct the agreement or disagreement of signs 
between all the particles. Finally we have slightly modified the construction presented in the last 
section by drawing from a symmetric Bernoulli the sign X{H) assigned to the H-ih. particle rather 
than the sign of a fixed particle: this clearly will not change the final distribution at all. All the 
elements of the desired vector of signs can then be obtained by multiplying Z L for a column vector 
with all zeros except X(H) in the H-th position. 

With a few changes, the algorithm can be adapted to the biased case in which the term ^ is 
added to the r.h.s. of ([2]). Recall that in this case ro = YlT=i \@k\ < 1 so that the mass 1 — ro is 
placed on to get a mass distribution on the integers 9* which replaces 9. Now at each step there 
is a positive probability that a particle "gets stuck" in a site. When the J-th particle stops its 
sign is sampled from the distribution of a binary random variable U with P(U = ±1) = ^z^y ■ 
This is sampled independently for particles stuck in different positions. The final vector of signs 
is obtained by post-multiplying Z L for the column vector containing the signs drawn for the stuck 
particles and zeros in the other positions. Such an algorithm is a slight modification, in order to 
profit of coalescence, to the one described in [CFF02 . 

However, in the latter paper, a general purpose "coupling from the past" algorithm was de- 
scribed, based on general bounds on the "dependence of the kernel on the past". For this general 
algorithm to work, in the case YlkLi \@k\ < 1, a control on the tails of the distribution 9* is needed, 
which is quite stronger than that required by Theorem El We should also mention that in [D SP12| 
a modification of such an algorithm is proposed, which works also for an example in which 9q = 
and EfcLi \9k\ = 1. But such a modification was rather ad hoc, and it was not clear how to pursue 
the same strategy for a generic vector 9 with YlkLi \@k\ = 1- 

It is worth to mention that in the case N > 1 a.s. it is still possible to produce a compatible 
process by repeating independently for each tree of the graph (Z, B) the random assignment of sign 
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to the two sets of the partition described above. But the uniqueness proof as well as the perfect 
simulation algorithm given here obviously fails, even if uniqueness can still hold. 

Finally, in dimension higher than one, the reader is referred to [G LO101 IDSL12] for perfect 
simulation of infinite range Gibbs fields. 
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